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Abstract 

', CLIFFORD performs various computations in GraBmann and Clifford al- 

gebras. It can compute with quaternions, octonions, and matrices with 
entries in Ci{B) - the Clifford algebra of a vector space V endowed with 
(Nl , an arbitrary bilinear form B. Two user-selectable algorithms for Clif- 

ford product are implemented: cmulNUM - based on Chevalley's recursive 
formula, and cmulRS - based on non-recursive Rota-Stein sausage. GraB- 
mann and Clifford bases can be used. Properties of reversion in undotted 
I and dotted wedge bases are discussed. 
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^ 1. Introduction 

I As many programs, CLIFFORD emerged from a practical problem. Relatively 

complicated algebraic manipulations with octonions, which can be performed in 
spin(7), started a thread which has now developed into a multi purpose algebra 
tool. It is the basic structure of a vector space V endowed with a quadratic 
form Q which is common to a vast host of mathematical, physical and engi- 
neering problems, and which, on the other hand, allows one to build naturally 
-i.e., in a categorial sense 'for free'- an algebra structure, the Clifford algebra 
Ci{V,Q). While in a conventional vector space framework one makes a good 
use of the vector space structure, it is tedious to compute with vectors since 
vector multiplication is lacking. Having established a Clifford algebra structure 
provides one with an entirely new formalism that now can be applied to solving 
completely different problems. 

In this sense, CLIFFORD is a basic tool for all such investigations and appli- 
cations which can be carried in finite dimensional vector spaces equipped with 

*The first author, R.A., acknowledges gratefully financial support from the College of Arts 
and Sciences, Tennessee Technological University. 
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a quadratic form or, equivalently, with a symmetric biUnear form commonly 
called inner or scalar product. The intrinsic abilities of Maple even allow to use 
CLIFFORD in projective and affine geometry, visualizing complicated incidence 
relations which are helpful, e.g., for image processing, visual perception and 
robotics. 

The authors of CLIFFORD have been interested in fundamental questions about 
g -deformed symmetries and quantum field theory. A reasonable number of new 
results has been derived by using systematically the ability to compute with a 
Computer Algebra System (CAS) at hand. Moreover, just asking questions as 
"what is the most general element fulfilling ..." has led to unexpected results 
and new insights. Testing of theorems, usually to check the consistency of the 
software, unfortunately has led from time to time to counter examples, that have 
made a rethinking and a reformulation necessary. However, the most striking 
ability of CLIFFORD is that it is unique in being able to handle Clifford algebras 
of arbitrary bilinear form not restricted by symmetry and not directly related 
to a quadratic form. It is now well known that such structures are related to 
Hopf algebraic twists, and the later versions of CLIFFORD make a good use of a 
process called Rota-Stein clijjordization, which turns out to be a Drinfeld twist 
of a GraBmann Hopf algebra ( [Brouder et al., 2002| ). 

The present paper introduces the reader to the package. It is assumed that 
she is already familiar with Maple ( |Maple, 2003|) , a general purpose CAS; if not 
please consult e.g. ([Wright, 2002|) . Of course, such a paper cannot be a user guide 
but may only be a demonstration of the usability and strength of the package. 
The interested reader is invited to download the package and to have a closer 
look at the online documentation in the Maple help browser or to download a 
pdf file of approx. 500 help pages. Therein we provide also further mathematical 
background and references and a detailed description for every function. How- 
ever, the present article provides also a quick -and dirty- introduction which is 
sufficient to get started. 

A list of aims behind experimental mathematics which are the guiding beacons 
for this work is found in a companion paper ( [Ablamowicz and Fauser, 2003b ) 
that describes the supplementary package BIGEBRA. 



2. Notations and basic computations 

CLIFFORD uses as default a standard GraBmann basis (GraBmann multivectors) 
in /\V where V = span {ej 1 1 < i < n} for 1 < n < 9. Then /\V = 
span {1, Bi A ej A . . . A ek\l < i,j,k < n}. In CLIFFORD these basis monomials 
are normally written as strings {Id, el, ... , e9, elwe2, elweB, . . . , elwe2we3, . . . , 
etc.} although they can be aliased to shorten input. Here elwe2 is a string that 
denotes ei A 62 and Id denotes the unity 1 in /\V. For example, when n = 3, 
GraBmann basis monomials are: 
> W=cbasis (3) ; 

W = [Id, el, e2, e3, elwe2, elweS, e2we3, elwe2we3] 
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but aliases can also be used to shorten input /output: 

> eval (makealiases (3) ) ; 

/, el2, e21, el3, e31, e23, e32, el23, el32, e213, e231, e312, e321 
In the above, eijk = eiwejwek is the wedge product of three 1-vectors: e^, e^, e^. 
Thus, the most general element of the GraBmann algebra /\V is a. GraBmann 
polynomial which is just a linear combination of GraBmann basis monomials with 
real coefficients. Notice that symbolic indices are allowed: 

> pi : =Id+4 . 5*ei-alpha*elwe2we3; 

pi := Id + 4.5 ei-a el23 
Reordering of GraBmann monomials can be explicitly accomplished by the pro- 
cedure reorder while CLIFFORD procedures ordinarily return their results in 
standard (reordered) basis. 

> p2 : =-e3we2wel-xO*Id+xl2*e2wel+a*ejwei ; reorder (p2) ; 

p2 := -e321 - xO Id + xl2 e21 + a ejwei 
el23 - xO Id - xl2 el2 - a eiwej 
The wedge product A is computed with the procedure wedge or its ampersand 
counterpart : 

> wedge (el , e2) , el &w e2;wedge(ea,eb,ec) ,ea &w eb &w ec;pl &w p2; 

el2, el2 
eawebwec, eawebwec 
el23 -xOId- 4.500000000 xO el + a xO el23 - xl2 el2 
The Clifford product can be introduced in /\V by means of a left J b (or 
right Lfi) contraction dependent on an arbitrary bilinear form B : V x V W 
following Chevalley's recursive definition and it will be discussed in Sect. ^ This 
leads to elements of the Clifford algebra C£{B) expanded into multivectors. 
Clifford multiplication is then implicitly dependent on B. It is given by the 
procedure cmul with the infix form &c. 

> cmul(el,e2) ,&c(el,e2) ;cmul(ea,eb,ec) ; 

el2 + Bi^2ld, el2 + Bi^2ld 

eawebwec + Bh^c e« — Ba^c ^b + Ba^b ec 
Simultaneous computation in C£{K) and C£{B) can be performed since cmul 
can accept a name of a bilinear form as a parameter. For example, 

> cmul [K] (el , e2) , &c [K] (el , e2) ; cmul [K] (ei , e j , ek) ; 

el2 + Ki^2ld,el2 + Ki^2ld 
eiwejwek + Kj^k ei — Ki^k ej + Kij ek 
Of course, the form B or K can be numeric or symbolic. For example, 

> B: =matrix(2,2, [1 ,a,a, 1] ) ; 

B:=\^ ?" 
a 1 

then GraBmann basis for Ci{B) ot /\V will be: 

> cbas : =cbasis (2) ; 

cbas := [Id, el, e2, el2] 
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while the Chfford multipHcation table of the basis GraBmann monomials will be 
as follows: 

> MultTable : =matrix(4,4, (i, j)->cniul(cbas [i] ,cbas [j] )) ; 



MultTahle : = 



Id 
el 

e2 
el2 



el 
Id 

-el2 + ald 
a el — e2 



e2 

el2 + ald 
Id 

el - ae2 



el2 
e2 - ael 
ae2 - el 
-l + a^)Id 



Of course, the Grafimann multiplication table will be: 

> wedgetable : =matrix(4,4, (i , j ) ->wedge (cbas [i] , cbas [j] ) ) ; 

Id el e2 el2 



wedgetable :- 



el el2 
e2 -el2 
el2 

Let B = g + F where g (F) is the symmetric (antisymmetric) part of B : 
> g,F:=matrix(2,2, [gll,gl2,gl2,g22] ) .matrix (2, 2, [0,F12,-F12,0] ) ; 
B : =evalm(g+F) ; 



gll gl2 
gl2 g22 





-F12 



F12 




B :-- 



gll gl2 + F12 
gl2 - F12 g22 
Then, the multiplication table of the basis monomials in Ci{B) will be: 

> MultTable : =matrix(4,4, (i , j )->cmul (cbas [i] , cbas [j] ) ) ; 

MultTable := 
[Id , el , e2, el2] 

[el , gll Id , el2 + {gl2 + F12) Id , gll e2 - {gl2 + F12) el] 
[e2 , {gl2 - F12) Id - el2 , g22 Id , {gl2 - F12) e2 - g22 el] 
[el2 , {gl2 - F12) el - gll e2 , g22 el - {gl2 + F12) e2 , 

{gl2'^ - F12^ - g22 gll ) Id - 2 el2 F12] 
Observe, that the "standard" relations e^ej + ejGj = {Bij + = '^Oij^ are 

satisfied by the generators ej, i = 1,2, ... ,n, irrespective of the presence of the 
antisymmetric part F in B. For example, 

> cmulEg] (el,e2)+cmul[g] (e2,el) ; 

2 Id gl2 

> cmul[B] (el,e2)+cmul[B] (e2,el) ; 

{gl2 + F12) Id + {gl2 - F12) Id 

> clisort (simplify (%)) ; 

2 gl2 Id 

It is well known [Lounesto (2001|) , [Lam (1973|) that real Clifford algebras 
Cl{y, Q) = Cip^q are classified in terms of the signature {p, q) of Q and the di- 
mension dim (V") = n = p + q. Information about all Clifford algebras Cip,q, 1 < 
n < 9 for any signature {p, q) has been pre-computed and stored in CLIFFORD. 
It can be retrieved with a procedure clidata. 
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> clidataC [2,0] ) ; #Clifford algebra of the Euclidean plane 

[Id, e2], [Id], [Id, e2]] 



[real, 2, simple, — Id + — el , 



The above output indicates that the CUfford algebra C£2 of the Euchdean plane 
is a simple algebra isomorphic to Mat(2,M) while the 4th entry in the data 
list is a primitive idempotent / then has been used to generate a minimal left 
spinor ideal S = Cl2.f and, subsequently, the left spinor (lowest dimensional and 
faithful) representation of CC.2- In general it is known that, depending on {p, q) 
and n = dim{V), the spinor ideal S = Cip^gf is a right /C-module where K 
is either M, C, or EI for simple Clifford algebras when {p — q) ^ I mod 4 or 
M © M and H © H for semisimple algebras when {p — q) = 1 mod 4 ( |Lounest"o 



198 Ij ), (Pelmstetter, 19821) . Elements in the third list [Id, e2] generate a real 
basis in S modulo /, that is, S = {IdSzc f,e2 Szc f) = (/, e2&c/). Elements 
in the 5th list generate a subalgebra F of Ci{Q) that is isomorphic to K. In 



the case of Ci2 we have that F = {Id) 



The last list gives 2 generators 



modulo f of S viewed as a right module over K where k = q — rg_p and r 
is the Radon-Hurwitz number.]] Number k gives a number of factors of type 
|l + |Tj, where {Ti}, i = 1, . . . , k, are commuting basis GraBmann monomials 
that square to 1 whose product gives a primitive idempotent / in Ci{Q). For 
all Clifford algebras C£{Q), l<n = p + q<9, and for any signature {p,q) 
of Q, spinor representations have been pre-computed ([Ablamowicz, 1998| ) and 
can be retrieved with the procedure matKrepr. For example, 1-vectors ei and 
62 in Ci2 are given in the (/, e2 &c /) basis of = €£2/ as: 
> matKrepr ( [2 , 0] ) ; 



[el 



e2 



In another example, Clifford algebra Ci^ of is isomorphic with Mat(2,C) : 
> B:=linalg[diag] (1,1,1) : clidata( [3,0] ) ; 



^Id + ^el, [Id, e2, e3 , e23\, [Id, e23], [Id, e2]] 



[complex, 2, simple 

and its spinor representation is given in terms of Pauli matrices: 
> matKrepr ( [3,0] ) ; 



[el 



e3 





e23 



-e23 




1 9 _ 1 
-1 J ' '''^ ~ [ 1 
Notice that F = {Id, e23) (e23 = e2we3) is a subalgebra of Cl^ isomorphic 
to C. Since Pauli matrices belong to Mat(2,F), it is necessary for CLIFFORD to 
compute with Clifford matrices, that is, matrices of type climatrix with entries 
in a Clifford algebra. 

> Ml , M2 , M3 : =rhs (% [1] ) , rhs (% [2] ) , rhs (% [3] ) ; 



Ml, M2, M3 



^Type TRHnumber in a Maple session when CLIFFORD is installed for more help. 



" 1 


■ 




"01' 




-e23 ' 





-1 




1 


•> 


e23 
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Of course Pauli matrices satisfy the same defining relations as the basis vectors 
61,62, and 63 :| For example: 

> 'Ml &cm M2 + M2 fecm Ml' = evalm(Ml &cm M2 + M2 fecm Ml); 
'el &c e2 + e2 &c el'=el &c e2 + e2 &c el; 



Ml kcm M2 + M2 kcm Ml 








el kc e2 + e2 kc el = 
> 'Ml &cm Ml' = evalmCMl &cm M1),'M2 &cm M2' = evalm(M2 &cm M2) , 
'M3 &cm M3' = evalm(M3 &cm MS); 

'el &c el' = el &c el,'e2 &c e2' = e2 &c e2,'e3 &c e3' = eS &c e3; 



Ml kcm Ml 



1 




M2 kcm M2 



M3 kcm M3 



el &c el = Id, e2 kc e2 = Id, eS kc eS = Id 
The procedure matKrepr gives the linear isomorphism between C£{Q) and 
Mat(2,M), and, in general, between Ci{Q) and Mat(2*^,K),K = M,C,H, for 
simple algebras and Ce{Q) and Mat(2^ K)®Mat{2'', K),K = R,m, for semisim- 
ple algebras. In this latter case, it is customary to represent an element in Ci{Q) 
in terms of a single matrix over a double field M © M or EI © HI rather than as 
pair of matrices.0 

One can easily list signatures of the quadratic form Q for which Ci{Q) is 
simple or semisimple. For more information, type ?all_sigs. For example, 
has a spinor representation given in terms of 2 by 2 quaternionic matrices whose 
entries belong to a subalgebra F of C£i,3 spanned by {Id, e2, e3, e2we3) : 

> B:=linalg[diag] (1,-1,-1,-1) : clidata( [1 ,3] ) ; 

[quaternionic, 2, simple, 2^^^^ 2 t^*^' ^'^^ ^-^^-^ ^^^5], 

[Id, e2, e3, e23], [Id, el]] 

> matKrepr ( [1,3] ) ; # quaternionic matrices 



[el 




1 



e2 



e2 






-e2 



e3 



e3 






-e3 



e4 



CLIFFORD has built-in several special-purpose procedures to deal with quater- 
nions and octonions (type ?quaternion and ?octonion for help). In particular, 
following [Lounesto et al. (1987| ), octonions are treated as para- vectors in C£oj 
while their non-associative multiplication, defined via Fano triples, is related to 
the Fano projective plane F2 (see ?omultable, ?Fano_triples for more infor- 
mation). Since the bilinear form B can be degenerate^], one can use CLIFFORD 
to perform computations in C^o,3,i, the Clifford algebra of the quadratic form 

■l-Here &cto is a matrix product where Clifford multiplications is applied to the matrix 
entries. See ?&cm for more information. 

''Procedures adf matrix and mdf matrix add and multiply matrices of type df matrix over 
such double fields. For more information see TmatKrepr. 

^When B = Q then C£(V, B) ~ /\V and computations in the Graf5mann algebra /\ V can 
then be done with CLIFFORD. 
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(5(x) = —xf — X2 — xl where x = xiei + 2:262 + x^e^ + x^e^^ G M^. This algebra 
is used in robotics to represent rigid motions in and screw motions in terms 
of dual quaternions ( |Selig, 19"96| ). 



Thus, CLIFFORD is a repository of mathematical knowledge about Clifford al- 
gebras of a quadratic form in dimensions 1 through 9. Together with the supple- 
mentary package BIGEBRA, ( [Ablamowicz and Fauser, 200 1|) it can be extended 



to graded tensor products of Clifford algebras in higher dimensions. Package 
BIGEBRA (see ?Bigebra) is described in a companion paper ([Ablamowicz and 



Fauser, 2003b| ). For more information about any CLIFFORD procedure, see its 
help page through the Maple browser. For a computation of spinor representa- 
tions with CLIFFORD, see [Ablamowicz (1998| ). 

3. Clifford product 

3.1. Recursive procedure cmulNUM 

Since the Clifford product provides the main functionality of the CLIFFORD pack- 
age, the best available mathematics has been used in its code. The user normally 
does not use the internal functions cmulRS and cmulNUM but the wrapper func- 
tion hc[K]{argl,arg2^ . . ) or, in long form, cmil[K]{argl,arg2, . . .) to pass 
the name of a bilinear form K. The wrapper function can also act on any num- 
ber of arguments of type clipolynom (since the Clifford product is associative 
this makes sense) and on a much wider class of types including Clifford matrices 
of type climatrix. It can also accept Clifford polynomials in other bases such 
as the Clifford basis {1, ej, ej &C e^-, &C(ej, e^, e^), . . .} where &C denotes the 
unevaluated Clifford product. Clifford basis differs from the Grafimann exterior 
basis when B is not a diagonal matrix.|| 

There are two internal Clifford multiplication procedures which are appropri- 
ate for different purposes: cmulRS and cmulNUM. While cmulNUM is fast on 
sparse numeric matrices and on numeric matrices in general for dimensions 
dim(1/) > 5, procedure cmulRS was designed for symbolic calculations. Since 
cmulRS computes reasonably well in the numeric sparse case up to dim {V) = 5, 
it was chosen as the default product of the package. Both internal Clifford mul- 
tiplication procedures take two Clifford monomials of type clibasmon as input 
together with a third argument of type name , symbol , matrix or array which 
represents the chosen bilinear form. 

There is a facility to let the user select cmulRS or cmulNUM when knowledge 
of the bilinear form allows one to decide which procedure might yield better 
performance. The user can supply a new product function (not necessarily a 
Clifford product) acting on two basis monomials. The wrapper uses an appropri- 
ate function, cmulRS, cmulNUM, or the other, which can be selected via a special 
function useproduct. 

''Procedures converting between Grafimann and Clifford bases are part of a supplementary 



package CliSplus (Ablamowicz and Fauser, 200f) while Clifford polynomials expressed in 



the Clifford basis are of type cliprod. Type ?cliprod for more information. 
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The evaluation of Clifford products in a GraBmann basis is quite involved and 
normally is done by a recursive process that involves Chevalley deformation. 
This approach has been employed in cmulNUM. However, Hopf algebraic methods 
can be apphed also, and have been used in cmulRS. Unfortunately, during the 
recursive evaluation many unnecessary terms are calculated which later cancel 
out at the next recursive call. This feature, while being beneficial when the 
bilinear form is sparse numeric since it cuts out many branches of the recursion 
quite early, prevents fast evaluation in the symbolic case where in general all 
terms might be non-zero. From this observation, the two possibilities to evaluate 
the Clifford product have emerged. 

We introduce the Chevalley deformation and the Clifford map to explain the 
algorithm used in cmulNUM. The Clifford map 7x is defined on u & /\V as 

(i) 7x(w) = LC{x, u, B) + wedge(x, u)=^-^BU + y./\u 

(ii) 7x7y = 7xAy + 5(x, y)7i 

(iii) 7ax+fey = a7x + &7y 

where x,y G 1/ ( [Lounesto, 200 1| ). One knows how to compute with the wedge 
X A -u and the left contraction LC(x, u, S) = x ?/with respect to the bilinear 
form B (in CLIFFORD left contraction is given by the procedure LC). Following 
Chevalley, the left contraction has the following properties: 

(i) xJBy = S(x,y) 

(ii) X ("U A t;) = (x m) A f + -u A (x f ) 

(iii) (u A f ) AbW = u-\b{v-^bw) 

where xgV^, ii,fG/\y and ii is GraBmann grade involution. Hence we can use 
the Clifford map 7x (Chevalley deformation of the GraBmann algebra) to define 
a Clifford product of a one- vector and a multivector. Analogous formulas can be 
given for a right Clifford map using the right contraction implemented as 
the procedure RC. 

The Clifford product cmul or its ampersand form &c can now be defined as 
follows: We have to split off recursively a single element from the first factor of 
the product of two GraBmann basis monomials and then use Chevalley's Clifford 
map. For example, 

(ea A . . . A eb A ej &c (e/ A . . . A e^) = 

(e„ A . . . A et) &c (e^ (e/ A . . . A e^) + A e/ A . . . A e^) 

- (e^ A . . . A Gfc ej &c (e^ A ... A e^). (1) 

Specifically, for (ei A 62) &c (es A 64) we have 

(ei A 62) &c (ea A 64) = (ei &c 62) &c (ea A 64) — B{ei, 62)! &c (ea A 64) 

= ei &c (5(62,63)64 - 5(62,64)63 + 62 A 63 A 64) 
-5(61,62)1 &c (63 A 64) 
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and a second recursion of the process gives now 

= B{e2, e3)B(ei, 64) - B{e2, 64)^(61, eg) + B{e2, e3)(ei A 64) 
-S(e2, e4)(ei A eg) + B(ei, e2)(e3 A e4) - B{ei, e3)(e2 A 64) 
+S(ei, e4)(e2 wes) + ei A e2 A 63 A 64 - B(ei, e2)(e3 A 64) 

with the bolded terms canceUing out. Note that the last term in the r.h.s. was 
superfluously generated in the first step of the recursion. 

The Clifford product can be derived from the above recursion by linearity and 
associativity. The induction starts with a left factor of grade one or grade zero 
which is trivial, i.e. 1 &c e„ A . . . A e^, = e^ A . . . A e^. In the case when the left 
factor is of grade one, we use the Clifford product expressed by the Clifford map 
of Chevalley, i.e., Bq &c e?, A . . . A ec = e^ e?, A . . . A e^ + e^ A ef, A . . . A ec. 
We make a complete induction in the following way: If the left factor is of higher 
grade, say n, one application of the recursion yields Clifford products where the 
left factor side is of grade either n — 1 or n — 2, hence the recursion stops after 
at most n — 1 steps. 

A disadvantage of the recursive approach is that additional terms are produced 
by shifting GraBmann wedge products into Clifford products to swap one factor 
to the right. These terms cancel out, but this process increases unnecessarily 
computing time. 

An advantage of the recursive approach occurs when the bilinear form B is 
numeric and sparse, that is, with many zeros. In this case after any recursion 
many terms drop out since Maple immediately simplifies such expressions and 
only a few remaining terms enter the next step of the recursion. If the dimension 
of V is large, i.e., dim (V) > 6, sparse matrices benefit drastically from this 
process over the Hopf algebraic approach of cmulRS which computes all terms 
without benefiting from the special sparse numeric form of the bilinear form. 

One could think about shifting factors from right to left, however this works 
out in the same way. Moreover, if the grade of the left factor n is greater than 
the grade m of the right factor, then the recursion stops also (since the terms 
evaluate to zero) after at most n — 1 steps, so no increase in performance can be 
gained this way. For example, the above computation in Maple will be performed 
as follows: 

> cmul(elwe2,e3we4) ; 

{B2, 3 Si, 4 - B2, 4 Si, 3) Id + B2, 3 el4 - B2, 4 el3 - Si, 3 e24 + Si, 4 e23 + el234 
Notice also that cmul accepts an arbitrary bilinear form K as its argument: 

> cmul [K] (elwe2,e3we4) ; 

{K2, 3 i^i, 4 - i^2, 4 i^i, 3) Id + K2, 3 eU - K2, 4 el3 - Ki, 3 e24 + i^i, 4 e23 + el234 
and likewise its ampersand formQ 

**Procedures cmulNUM and cmulRS do not have their special ampersand forms. However, 
procedure &c uses internally one procedure or the other depending on the current value of an 
environmental variable _def ault_Clif f ord_product. Current values of these variables can be 
displayed by CLIFFORD_ENV. 
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> &c [K] (ei , ejwekwel) ; 

eiwejwekwel + Kij ekwel — Ki^ ^ ejwel + i ejwek 
where we have also shown the abihty to use symbohc indices. For clarity and to 
show our approach we display the algorithm of cmulNUM in Appendix A. 



3.2. Procedure cmulRS based on the Rota-Stein combinatorial process 

The procedure cmulRS is computed using the non-recursive Rota-Stein clif- 



fordization ( [Rota and Stein, 1994| ; [I'auser, 2002| ; [Ablamowicz and l:''auser, 2003a| JED 
and the help pages of the BIGEBRA package for further literature. The cliffordiza- 
tion process is based on Hopf algebra theory. The Clifford product is obtained 
from the GraBmann wedge product and its GraBmann co-product as follows (in 
tangle notation): 




(2) 



where A is the GraBmann exterior wedge product and Aa is the GraBmann 
exterior co-product, which is obtained from the wedge product by categorial 
duality, i.e. to every algebra over a linear space A having a product we find a 
co-algebra having a co-product over the same space by reversing all arrows in 
all axiomatic commutative diagrams. Note that the co-product splits each input 
'factor' X into a sum of tensor products of ordered pairs X(i)j,X(2)j. The main 
requirement is that every such pair multiplies back to the input x when the 
dual operation of multiplication is applied, i.e. A X(2)i = x for each ith 
pair. The 'cup' like part of the tangle decorated with is the bilinear form 
B on the generating space V extended to the whole GraBmann algebra, i.e. the 
map B^ : /\V X /\V k with B : V x V ^ k as -B(x, y) on vectors. Hence, 
cmulRS computes on GraBmann basis monomials x and y for the given B , later 
extended to polynomials by bilinearity, as follows: 



cmulRS(x, y,B) = ^ ^(±)x(i), A y{2)jB{x(2)i,y{i)j) 
i=i j=i 



(3) 



where n, m are dummies describing the cardinalities of the required splits and 
the sign is due to the parity of a permutation needed to arrange the factors. 
The simplified algorithm of cmulRS is as follows: 



1 > cmulRS (x,y,B) [x, y two Grassmcinn monomials, B 

2 > begin 

3 > Istx <- list of indices from x 

4 > Isty <- list of indices from y 



bilinear form] 
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5 > NX <- length of Istx 

6 > NY <- length of Isty 

7 > funx <- function maps integers 1..NX onto elements of Istx keeping their order 

8 > funy <- function maps integers 1..NY onto elements of Isty keeping their order 

9 > (this is to calculate with arbitrary indices and to compute necessary signs) 

10 > psetx <- power set of 1..NX (actually a list in a certain order) 

11 > (the i-th eind (2~NX+l-i) -th element are disjoint adding up to the set 1..NX 

12 > psety <- power set of 1..NY (actually a list in a certain order) 

13 > (the i-th and (2"NY+l-i)-th element are disjoint adding up to the set 1..NY 

14 > (for faster computation we sort this power sets by grade) 
16 > (we compute the sign for any term in the power set) 

16 > psetx <- sort psetx by grade 

17 > psety <- sort psety by grade 

18 > pSgnx <- sum_(i in psetx) (-l)~sum_(j in psetx[i]) (psetx[i] [j]-j) 

19 > pSgny <- sum_(i in psety) (-l)"sum_(j in psety[i]) (psety [i] [j] -j ) 

20 > (we need a subroutine for cup tangle computing the bilinear form cup(x,y,B)) 



21 > begin 

22 > if |x| <> lyl then RETURN (0) fi 

23 > if Ixl = then RETURN(l) fi 

24 > if 1x1 = 1 then RETURN (B [x [1] ,y [1] ] ) fi 

26 > RETURN(sum_(j in 1 . . |x| ) (-1) " (j-l)*B(x[l] ,y [j] )*cup(x[2. .-1] ,y/y [j] ,B) ) 

26 > end cup 



27 > (now we compute the double sum, to gain efficiency we do this grade wise) 

28 > (note that there are r over NX r-vectors in psetx, einalogously for psety) 

29 > max_grade <- I Istx <- convert_to_set union Isty <- convert_to_set I 

30 > res <- 0, posl <- 

31 > for j from to NX (iterate over all j -vectors of psetx) 



32 > begin 

33 > FI <- Nl!/((Nl-j) !*j !) (number of terms (Nl over j)) 

34 > pos2 <- 

35 > for i from to min(N2 ,max_grade-j ) 

36 > (iterate over all i-vectors of psety not exceeding max_grade while 

37 > others are zero) 
,58 > begin 

39 > F2 <- N2!/((N2-i) !*i!) (number of terms (N2 over i)) 

40 > for n from 1 to FI (for all j-vectors) 

41 > begin 

42 > for m from 1 to F2 (for all i-vectors) 

43 > begin 

44 > res <- res + 

46 > pSgnx [posl+n] *pSgny [pos2+m] * 

46 > *cup(funl(psetx[PNl-posl-n] ) ,fun2 (psety [pos2+m] ) ,lname)* 

47 > makeclibasmon -> ( [f unl -> psetx [posl+n] ,fun2 -> psety [PN2-pos2-m] )] ) 

48 > end 

49 > end 

50 > pos2 <- pos2+F2 

61 > end 

62 > posl <- posl+Fl; 

53 > end 

64 > reorder -> res (reorder basis elements in res into standard order) 



65 > end cmulRS 



It is clear from this algorithm that only those terms are considered which might 
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be non-zero: If all Bij are non-zero and different so that no cancellation takes 
place between them, all these terms will survive. The combinatorial power of 
the Hopf algebraic approach is clearly demonstrated with this algorithm and its 
superior behavior shows up in benchmarks ( [Ablamowicz and Fauser, 2003a ). 



4. Dotted and undotted Grafimann bases in quantum Clif- 
ford algebras 

It was shown above that CLIFFORD uses the GraBmann algebra /\V as the 
underlying vector space of the Clifford algebra C£{V,B). Thus, the GraBmann 
wedge basis of monomials is the standard basis used in CLIFFORD. A general 
element u in C£{V, B) can be therefore viewed as a GraBmann polynomial. 

When the bilinear form B has an antisymmetric part F = —F'^, it is conve- 
nient to split B = g + F, where g is the symmetric part of B, and to introduce 
the so called "dotted GraBmann basis" ( [Ablamowicz and Lounesto, 1996]) and 



the dotted wedge product A. The old GraBmann basis will be referred to here 
as "undotted GraBmann basis". In CLIFFORD, the wedge product is given by the 
procedure wedge and while the dotted wedge product is given by dwedge 
and Szdw. 

According to the Chevalley definition of the Clifford product &c, we have 

x&c 'U = xJbm + x&w u = LC(x, u, B) + wedge (x, u) (4) 

for a 1-vector x and an arbitrary element u of Ci{B). Here, LC(x, -u, S) de- 
notes the left contraction of m by x with respect to the entire bilinear form B. 
However, when B = g + F then 

X M = LC(x, u,B) = yi-igU + yi-iFU = LC(x, u, g) + LC(x, w, F) (5) 

and 

X &c M = LC(x, -u, S) + X &w u (6) 

= LC(x, M,5() + LC(x, u,F) + x &w u (7) 

= LC(x, M,5() + x &dw M (8) 

where x &dw m — x &w u = LC(x, F). That is, the wedge and the dotted wedge 
"differ" by the contraction term(s) with respect to the antisymmetric part F 
of B. This dotted wedge ^dw can be extended to elements of higher grades. 
Its properties are discussed next. 



4.1. Indexing dwedge and hdw 

Procedure dwedge (and its infix form ^dw) requires an index which can be a 
symbol or an antisymmetric matrix. That is, dwedge computes the dotted wedge 
product of two GraBmann polynomials and expresses its answer in the undotted 



Rafal Ablamowicz and Bertfried Fauser: Mathematics of CLIFFORD 



13 



basis. Special procedures exist which convert polynomials between the undotted 
and dotted bases. When no index is used, the default is F : 

> dwedge[K] (el+2*e2we3,e4+3*elwe2) ;&dw(ei+2*ejwek,ei+2*ejwek) ; 

-(-i^i,4 + 6 i^2, 3^^1,2) /rf - 6Ki^2e2we3 - 6X2,3 elwe2 - 2X2,4 e3 + 2X^,4 62 
- 3i^i,2 ei + elwe4 + 2 e2we3we4 

4 eiwejwek — 4Fi,k ej + 4 Fi, j ek — 8 Fj, k ejwek — 4 Fj, Id 
Observe that conversion from the undotted wedge basis to the dotted wedge 
basis using antisymmetric form F and dwedge[F] are related through the fol- 
lowing convert function: 

dwedge[F](el, e2, ...,en) = convert(eliye2w...tuen, wedge_to_dwedge, F) 

which can be shown as follows: 

> F : =array (1 . .9,1. . 9 , antisymmetric) : 

> dwedge [F] (el , e2)=convert (wedge(el ,e2) , wedge_to_dwedge ,F) ; 

elwe2 + Fi, 2 Id = elwe2 + Fi, 2 Id 

> dwedge [F] (el , e2 , e3)=convert (wedge (el , e2, eS) , wedge_to_dwedge,F) ; 

elwe2we3 + ^2,3 ei - Fi,3 e2 + Fi,2 e3 = elwe2we3 + ^2,3 el - Fi,3 e2 + Fi,2 e3 

4.2. Associativity of dwedge 

Operation dwedge is associative with the unity 1 = Id as a unit: 

> evalb (dwedge [F] (dwedge [F] (el,e2) ,e3)=dwedge [F] (el , dwedge [F] (e2,e3))) ; 

true 

For some arbitrary random Clifford polynomials^ u,v, z expressed in GraBmann 
undotted basis we can show associativity as follows: 

> u : =2*Id+el-3*e2we3 : v : =3*Id-4*elwe3+e7 : z : =4*Id-2*e3+elwe2we3 : 

> evalb (dwedge [F] (Id,u)=u) , evalb (dwedge [F] (u, Id)=u) ; 

true, true 

> evalb (dwedge [F] (dwedge [F] (u , v) , z) =dwedge [F] (u , dwedge [F] (v , z) ) ) ; 

true 

We have, therefore, the following identity that is satisfied by any two elements 
u and V in Ci{B), B = g + F, that are, by default, expressed in terms of the 
undotted GraBmann basis: 

u A V = {up A Vf)-f- (9) 

Here up and vp are the elements u and v expressed in the dotted basis with 
respect to the form F while (. . .)-f denotes conversion back from the dotted 
basis to the undotted basis w.r.t. —F = F^ . This can be illustrated in CLIFFORD 
as follows: 

t'l'In CLIFFORD ver. 6 and higher there are three procedures useful for testing that return 
a random GraBmann basis monomial, a random monomial and a random polynomial, respec- 
tively. See ?rd_clibasmon, ?rd_climon, ?rd_clipolynom. 
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> uu:=convert(u,wedge_to_dwedge,F) ; #u converted to dotted basis 
vv:=convert(v,wedge_to_dwedge,F) ; #v converted to dotted basis 

uu :=el - 3 e2we3 - 8^2,3/^ + 2Id 

vv := 3Id -4elwe3 -4Fi^3ld+ e7 

> outl :=dwedge [F] (uu, vv) : #dwedge computed w . r . t . F 

> out2:=convert(outl ,dwedge_to_wedge,-F) ; #back to undotted basis 

out2 := 3 ei - 9 e2we3 + 6 Id -8 elweS + elwe? - 3 e2we3we7 + 2e7 

> out3:=wedge(u,v) ; #direct computation of wedge product 

outs := 3 ei - 9 e2we3 + 6 Id -8 elwe3 + elwe? - 3 e2we3we7 + 2e7 
and it can be seen that out2 = out3. 

4.3. Dotted and undotted wedge bases 

The default GraBmann basis in Ci{B) used in CLIFFORD is undotted. How- 
ever, one can easily use the dotted basis. For example, we expand the basis of 
the original wedge into the dotted wedge and back. For this purpose we choose 
dim {V) = 3 and consider C£{B) with the antisymmetric part F. The undotted 
wedge basis for /\V is then: 

> w_bas :=cbasis(dim_V) ; ## the wedge basis 

w.bas := [Id, el, e2, e3, elwe2, elwe3, e2we3, elwe2we3] 

Now we map the convert function onto this basis to get the dotted wedge basis: 

> d_bas:=map(convert,w_bas,wedge_to_dwedge,F) ; 
test_wbas : =map (convert , d_bas , dwedge_to_wedge , -F) ; 

d.bas := [Id, el, e2, e3, elwe2 + Fi^2 Id, elwe3 + Fi^^Id, e2we3 + F2,zld, 

elwe2we3 + ^2,3 el - Fi,3 + Fi,2 e3] 

test-wbas := [Id, el, e2, e3, elwe2, elwe3, e2we3, elwe2we3\ 
Notice that only the unity 1 and the one vector basis elements remain unal- 
tered and that the other basis elements of higher grades pick up additional terms 
of lower grades (which preserves the filtration). It is possible to define aliases 
in CLIFFORD for the dotted wedge basis "monomials" similar to the GraBmann 
basis monomials. For example, we could denote the element elwe2 -\- F[l, 2] * Id 
by elWe2 (= ei A 62) and similarly for other elements: 

> alias (elWe2=elwe2 + F [1 ,2] *Id,elWe3=elwe3 + F[l,3]*Id, 
e2We3=e2we3 + F[2,3]*Id, 

elWe2We3=elwe2we3+F [2 , 3] *el-F [1 , 3] *e2+F [1 , 2] *e3) ; 

/, elWe2, elWeS, e2We3, elWe2We3 
and then Maple will display automatically dotted basis in terms of the aliases: 

> d_bas ; 

[Id, el, e2, e3, elWe2, elWe3, e2We3, elWe2We3] 
That is, as linear spaces we find isomorphisms 

Ci{B) ^ (1,61,62,63,61 A e2,ei A 63,62 A eg, ei A 62 A 63) 
= (1, 61, 62, 63, 61 A 62, 61 A 63, 62 A 63, 61 A 62 A 63) 

where 61 A 62 = elWe2, etc. 
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4.4. Contraction in dotted and undotted bases 

The contraction -ix w.r.t. any bilinear form K works on both undotted and 
dotted bases in the same manner which can be seen if we re-convert the dotted 
basis after the computation into the wedge (undotted) basis. In a reasonable 
setting, the antisymmetric bilinear form F would be the antisymmetric part 
of B. To read more about the left contraction LC in Ci{B) check the help page 
for LC or see ([Ablamowicz and Lounesto, 1996| ). We have the following identity 
for any two elements u and v in C£{B) expressed in the undotted GraBmann 
basis: 

V -iBU= {v -is uf)-f (10) 

As before, up is the element u expressed in the dotted basis while (. . .)_f 
accomplishes conversion back to the undotted basis. To illustrate this fact, we 
first contract from the left an arbitrary element u in C£{B) by 1, e,, A e^, Gj A 
Gj A efc (here we limit our example to dim iV) = < i,j,k < 3) and then we 
extend it to a left contraction by an arbitrary element v in Ci{B). 

> u: =add(x . i*w_bas [i+1] , i=0 . . 7) :uF: =convert (uw,wedge_to_dwedge ,F) : 
v:=add(y. i*w_bas [i+1] ,i=0. .7) : 

Contraction with respect to 1 : 

> evalb(LC(Id,u,B)=convert(LC(Id,uF,B) ,dwedge_to_wedge , -F) ) ; 

true 

Contraction with respect to e-j : 

> evalb(LC(ei,u,B)=convert(LC(ei,uF,B) ,dwedge_to_wedge , -F) ) ; 

true 

Contraction with respect to A Sj : 

> evalb(LC(eiwej ,u,B)=convert (LC(eiwej ,uF,B) ,duedge_to_wedge ,-F) ) ; 

true 

Contraction with respect to ej A ej A : 

> evalb(LC(eiwejwek,u,B)=coiivert(LC(eiwejwek,uF,B) ,dwedge_to_wedge,-F)) ; 

true 

Finally, contraction with respect to an arbitrary element v : 

> evalb (LC (vw , uw , B) =convert (LC (vw , uF , B) , dwedge_to_wedge , -F) ) ; 

true 



4.5. Clifford product in dotted and undotted bases 

We can build a Clifford algebra Cl{B) over each basis set, that is, over the undotted 
or dotted GraBmann basis, but with different bilinear forms: B = g or B = g + F re- 
spectively (following notation from ( [Ablamowicz and Lounesto, 1996 )). Let's compute 
various Clifford products with respect to the symmetric form g and with respect to 
the full form B using procedure cmul that takes a bilinear form as its index. As an 
example, we will use two most general elements u and v m. f\V when dim iy) = 3. 
Most output will be eliminated. 

> u:=add(x.k*w_bas [k+1] ,k=0. .7) :v:=add(y .k*u_bas [k+1] ,k=0. .7) : 
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We can then define in f\V Clifford product cmul [g] with respect to the symmetric 
part g and another Chfford product cmul [B] with respect to the entire form B : 

> cmulg:=proc() RETURN ( cmul [g] (args) ) end: 
cmulB:=proc() RETURN ( cmul [B] (args) ) end: 

Thus, we are ready to perform computations around our next commutative diagram, 
however output will be eliminated to save space. First, we compute Clifford product 
cmul [g] {u, v) in Ci{g) in undotted Grafimann basis. 

> uv:=cmulg(u,v) : ## Clifford product w.r.t. g in Cl(g) in wedge basis 

Now, we convert u and v to up and vp, respectively, expressed in the dotted wedge 
basis: 

> uF:=convert(u,wedge_to_dwedge,F) :vF:=convert(v,wedge_to_dwedge,F) : 

We now compute the Clifford product of up and fi? in Ci{B) in the dotted wedge 
basis, 

> uFvF:=cmulB(uF,vF) : ## Clifford product in C1(B) in dwedge basis 

convert back the above result back to the undotted wedge basis: 

> uv2 : =convert (uFvF,dwedge_to_wedge ,-F) : ## convert result dwedge->wedge 
and verify that the results are the same: 

> simplify (uv-uv2) ; ## show equality ! 



Thus, we have shown that the following identity involving cmul [g] and cmul [B] is 
true (at least when dim(y) = 3) :0 

{uv)g = uhcg V = (up &:cB vp)-p = {{up vp)b)~f (H) 

This shows that the Clifford algebra Ci{g) of the symmetric part g of B using the 
undotted exterior basis is isomorphic, as an associative algebra, to the Clifford algebra 
Ci{B) of the entire bilinear form B = g -\- F spanned by the dotted wedge basis if 
the antisymmetric part F of i? is exactly the same as F is used to connect the two 
basis. 

4.6. Reversion in dotted and undotted bases 

We proceed to show that the expansion of the Clifford basis elements into the dotted 
or undotted exterior products has also implications for other well known operations 
such as the Clifford reversion anti-automorphism ~ : Cl{B) C£{B), uv i— > vu, which 
preserves the grades in /\V [but not in /\V unless B is symmetric] Only when the 
bilinear form is symmetric, we find that the reversion is grade preserving, otherwise it 
reflects only the filtration: That is, reversed elements are in general sums of terms of 
the same and lower degrees. 

> reversion(elwe2,B) ; #reversion with respect to B 
reversion(elwe2 ,g) ; #reversion with respect to g (classical result) 

-elwe2 - 2Fi,2/(i 

— elwe2 

^Here, {uv)g is the Clifford product with respect to g while up^csvp and {ufVf)b are 
the Clifford products with respect to B, that is, in C£{g) and C£{B), respectively. 
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Observe in the above that only when Bi^2 = B2^i, the result — ei A e2 known from 
the theory of classical Clifford algebras is obtained. Likewise, 

> cbas : =cbasis (3) ; 

cbas := [Id, el, e2, eS, elwe2, elweS, e2we3 , elwe2we3] 

> map(reversion,cbas,B) ; 

[Id, el, e2, e3 , -elwe2 - 2Fi,2 -elweS -2Fi^^Id, -e2we3 - 2^2,3 /ti, 

-2^2,3 ei +2Fi,3e5 -2Fi,2e5 - elwe2we3] 
If instead of B we use a symmetric matrix g = (or the symmetric part of B), 
then 

> map (reversion , cbas , g) ; 

[Id, el, e2, e3, -elwe2, -elwe3, -e2we3, -elwe2we3] 
Convert now ei A 62 to the dotted basis to get ei A e2 = elWe2 : 

> convert (elwe2,wedge_to_dwedge,F) ; 

elWe2 

Apply reversion to elWe2 with respect to F to get the reversed element in the dotted 
basis: 

> reversed_elWe2:=reversion(elWe2,F) ; 

reversed-elWe2 := —elwe2 — Fi^2ld 
Observe, that the above element is equal to the negative of elWe2 just like reversing 
elwe2 with respect to the symmetric part g oi B : 

> reversed_elWe2+elWe2; 



Finally, convert reversed elWe2 to the undotted standard Grafimann basis to get 
—elwe2 : 

> convert (reversed_elWe2 ,dwedge_to_wedge , -F) ; 

— elwe2 

The above, of course, can be obtained by applying reversion to elwe2 with respect to 
the symmetric part g of B : 

> reversion(elwe2 ,g) ; #reversion w.r.t. the symmetric part g 

— elwe2 

This shows that the dotted wedge basis is the particular basis which is stable under the 
Clifford reversion computed with respect to F, the antisymmetric part of the bilinear 
form B. This requirement allows one to distinguish Clifford algebras C£{g) which have 
a symmetric bilinear form g from those which do not have such symmetric bilinear 
form but a more general form B instead. We call the former classical Clifford 
algebras while we use the term quantum ClifTord algebras for the general not 
necessarily symmetric case ( [Ablamowicz and Fauser, 2000 ). 



5. Conclusions 

This paper continues with the second part (Ablamowicz and Fauser, 2003b) about 
BIGEBRA where further aims and outlooks for the future applications of CLIFFORD 
and BIGEBRA are given. 
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Appendix A: Code of cmulNUM 

Here is a shortened code of the recursive procedure cmulNUM. 



56 > cmulNUM ( al ,a2,B) [al , a2 - two Grassmsinn monomials, B - name of bilinear form] 

57 > begin 

58 > if nargs<>3 then ERROR (' exactly three arguments are needed') fi: 

59 > if has(0,map(simplify, [al,a2])) then RETURN(O) fi; 

60 > if a2='Id' then RETURN(al) fi: 

61 > if al='Id' then RETURN(a2) fi: 

62 > L <- indices from al 

63 > N <- length of L 

64 > coB,ncmieB <- coefficient of B, B [to handle -B] 

65 > if N=0 then RETURN (coeff(al, Id) *a2) elif N=l then 

66 > L2 <- list of indices from a2 

67 > RETURN (reorder] (simplify (makeclibasmonC [L [1] ,op(L2)]) 

68 > +add((-l)~(i-l)*coB*nameB[L[l] ,L2[i]]* 

69 > makeclibasmon(subs(L2[i]=NULL,L2)) ,i=l. .nops(L2))))) ; 

70 > elif N=2 then 

71 > xl <- substring (al , 1 .. 2) : 

72 > x2 <- substring(al ,4 . . 5) ; 

73 > p2 <- procname(x2,a2,B) : 

74 > S <- clibilinear(xl,p2, procname,B) ; 

75 > RETURN(simplify(S-coB*nameB[op(L)]*a2)) 

76 > fi; 

77 > X <-cat(e,L[-l] ) ; 

78 > pl<-substring(al , 1 . . (3*N-4) ) ; 

79 > p2<- procname(x,a2,B) : 

80 > S<-clibilinear (pi ,p2, procname.B) 

81 > -add((-l)"(i)*coB*nameB[L[-i] ,L[-1]]* 

82 > procname(makeclibasmon(subs (L [-i] =NULL,L [1 . . -2] ) ) , a2,B) , i=2 . . N) ; 

83 > RETURNreorder(simplify(S))) ; 

84 > end cmulNUM: 
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